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ABSTRACT 

We perform a combined analysis of two stringent constraints on two Higgs doublet 
model, coming from the recently announced CLEO II bound on B(6 -> 57) and from 
the recent LEP data on the ratio r(z^hadrlns) • ^ e include one- loop vertex corrections 
to Z — > bb in the model. We find that although the CLEO II bound serves as the 
strongest constraint present in the charged Higgs sector of the model, the current 
LEP value for Rb may also provide a further constraint for tan/? < 1. 
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Despite the remarkable successes of the Standard Model(SM) in its complete 
agreement with current all experimental data, there is still no experimental informa- 
tion on the nature of its Higgs sector. The 2 Higgs doublet model(2HDM) is one 
of the mildest extensions of the SM, which has been consistent with experimental 
data. In the 2HDM to be considered here, the Higgs sector consists of 2 doublets, <pi 
and 02, coupled to the charge -1/3 and +2/3 quarks, respectively, which will ensure 
the absence of Flavor- Changing Yukawa couplings at the tree level [||]. The physical 
Higgs spectrum of the model includes two CP-even neutral Higgs(if°, h°), one CP- 
odd neutral Higgs(A°) , and a pair of charged Higgs(if ± ). In addition to the masses 
of these Higgs, there is another free parameter in the model, which is tan/? = t^/fi, 
the ratio of the vacuum expectation values of both doublets. 

With a revived interest on the flavor-changing-neutral-current (FCNC) b — > 57 
decay, spurred by the CLEO bound B{b -» sj) < 8.4 x 10" 4 at 90% C.L. [|, it was 
pointed out recently that the CLEO bound can be violated due to the charged Higgs 
contribution in the 2HDM and the Minimal Supersymmetric Standard Model(MSSM) 
basically if m H ± is too light, excluding large portion of the charged Higgs parameter 
space [|]. The recently announced CLEO II bound B(b — > 57) < 5.4 x 10" 4 at 95%0 
excludes even larger portion of the parameter space 0. It has certainly proven that 
this particular decay mode can provide more stringent constraint on new physics 
beyond SM than any other experiments]^. In this report, we will show that in 
addition to the CLEO II bound on b — > 57, one may be able to constrain the 2HDM 
further, by incorporating one-loop vertex corrections to Z — > bb, with the recent LEP 
data on which is defined to be the ratio Y{z^rmdrins) • ^ number of authors in 

fact have studied the one- loop corrections to T(Z — > bb) in the context of the SM|J 
and a few extensions of it [110, [IT], O although the emphases were different. 

In the 2HDM, b —> 57 decay receives contributions from penguin diagrams 
with W — t loop and — t loop. The expression used for B{b — > 57) is given by 



B{b - s 7 ) 6a [v W/23 A, + |(^ 4 /23 _ ^16/23)^ + c 
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where 77 = a s (M z ) / a s {m b ), I is the phase-space factor I(x) = 1 — 8x 2 + 8x 6 — 
x 8 — 24x 4 lnx, and f(m c /m b ) = 2.41 the QCD correction factor for the semileptonic 
decay. We use the 3-loop expressions for a s and choose Aqcd to obtain a s [Mz) 
consistent with the recent measurements at LEP. In our computations we have used: 
a s (M z ) = 0.118, B{b -> cev) = 10.7%, m b = 4.8 GeV, and m c /m b = 0.3. The Ay, A g 
are the coefficients of the effective bs^ and bsg penguin operators evaluated at the 
scale Mz- The contributions to A 1 ^ g from the W ± — t loop, the H ± — t loop are given 
in Ref|TB[. As mentioned above, the CLEO II bound excludes a large portion of the 



parameter space. In Fig. 1 we present the excluded regions in (m H ± } tan/3)-plane 
for m t = 120, 130, and 150 GeV, which lie to the left of each curve. We have imposed 
in the figure also the lower bound on tan/5 from ^ ^ tan/5 ^ — obtained by 
demanding that the theory remain perturbative[14]. We see from the figure that at 
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large tan /3 one can obtain a lower bound on m#± for each value of m^. And we obtain 
the bounds , m H ± > 160, 186, 244 GeV for m t = 120, 130, 150 GeV, respectively. Our 
strategy is now to sample tan/3 and m H ± from the allowed regions in Fig. 1, which 
will be in turn used below to calculate Rb at one-loop level in this model. 

In the SM, the diagrams for the vertex corrections to Z — ► bb involve top quarks 
and W ± bosons. However, in the 2HDM there are additional diagrams involving if ± 
bosons instead of W ± bosons. These additional diagrams have been calculated in 



Ref |10|, [11], 0] • Here we use the formulas given in Ref |TT[ . The calculation involves 



numerical evaluation of the reduced Passarino-Veltman functions! 151. Our numerical 



results agree with those in Fig. 4-5 in Ref|TT[. In our calculation, we neglect the 
neutral Higgs contributions which are all proportional to ml tan /3 2 and become sizable 
only for tan (3 > ^ and very light neutral Higgs ^ 50 GeV, but decreases rapidly 



to get negligibly small as the Higgs masses become ^ 100GeV|12|. We also neglect 



oblique corrections from the Higgs bosons just to avoid introducing more paramters. 
This correction grows as m 2 H ± if m#± 3> iriH ,h ,A , an d gives only a —0.1% correction 
for m H ± = 500 GeV. Although tan (5 1 seems more appealing because of apparent 
hierarchy m t 3> nib, there are still no convincing arguments against tan/3 < 1. Here 
we choose to explore the region of tan f3 ^ 1. In Fig. 2 we show the model predictions 
for Rj, as a function of m t in comparison with the SM prediction. The parameters 
were sampled from the allowed regions in the Fig. 1. Two horizontal solid lines in the 
figure represent the lower limits from the recent LEP data Rb = 0.2203 ± 0.00270. 
The upper one corresponds to the 1-cr value, and the lower one to the 1.64-cr value. 
We note that the deviations from the SM can be quite large for tan/3 < 1 because 
the charged Higgs contribution grows as m\j tan 2 /3 for tan/3 <C — . The deviations 
for tan/3 < 1 can be as large as —2.2% for m t = 150 GeV while they become much 
smaller for tan/3 > 1[[RJ. We have also considered other constraints from low-energy 
data primarily in I? — B,D — D,K — K mixing that exclude low values of tan /3|H], [IT]] . 
But tan (3 ^ 0.5 seems to be still allowed by these constraints for m t ^ 150 GeV and 
rriH± ^ 250 GeV. However, for small values of tan /3, as is seen from Fig. 2, the model 
predicts m t ^ 110 GeV at l-a level, which is in conflict with the recent bound from 
CDF||18|, m t > 113 GeV although the corresponding prediction at 1.64-a level is still 
allowed . These low values of tan/3 ^ 0.5 might be at the verge of being disfavored 
even at 90%C.L. as the experimental lower bound on m t tends to grow. Nevertheless, 
the CLEO II bound is still by far the strongest constraint present in the charged 
Higgs sector of the model especially for tan/3 <; 1. 

In the MSSM, the situation becomes much more complicated. In the calcula- 
tion of B(b — > 57) in Ref|4[], other important contributions such as chargino-squark 
contribution were not included. However, it was shown very recently [13|] that the 
bound on m H ± in Ref 0] could well be evaded in the full supersymmetric calculation 
since the b — >■ S7 amplitude vanishes in the exact supersymmetric limit. Therefore, 
with full calculations, one may get significantly smaller bounds on m#± than the 
ones given above |7J. However, adding supersymmetric particles raises Rb significantly 
above the SM value due to chargino-squark loops and neutralino-squark loops can- 
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celling out to a great extent the charged Higgs constribution and also the standard 



contribution |TT[]. Thus, this makes it extremely difficult for one to be able to constrain 
the MSSM with the current LEP data. 

In conclusion, we have performed a combined analysis of two stringent con- 
straints on the 2 Higgs doublet model, coming from the recently announced CLEO II 
bound on B(b —> 57) and from the recent LEP data on Rb. We have included one-loop 
vertex corrections to Z — > bb in the model. We find that although the CLEO II bound 
serves as the strongest constraint present in the charged Higgs sector of the model, 
the current LEP value for Rb may also provide a further constraint for tan/3 < 1. 
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Figure Captions 

Figure 1: The regions in (m H ±, tan (3) plane excluded by the CLEO II bound B{b — > 
S7) < 5.4 x 10~ 4 at 95%, for m t = 120, 130, 150 GeV in 2HDM. The excluded 
regions lie to the left of each curve. The values of m t used are as indicated. 

Figure 2: R b as a function of top mass in the SM(Solid), and the 2HDM for five 
different sets of tan (3 and m H ± sampled from the allowed regions in the Fig. 
1. The values of (tan/5, m H ±) used £1X6 clS indicated under each curve. Two 
horizontal solid lines represent the lower limits from the recent LEP data R b = 
0.2203 ±0.0027. The upper one corresponds to the 1-a value, and the lower one 
to the 1.64-cr value. Anything above these lines is allowed. 
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